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Abstract — We define a notion of network capacity region of 
networlts tliat generalizes tlie notion of networlt capacity defined 
by Cannons et al. and prove its notable properties such as 
closedness, boundedness and convexity when the finite field is 
fixed. We show that the network routing capacity region is a 
computable rational polytope and provide exact algorithms and 
approximation heuristics for computing the region. We define 
the semi-network linear coding capacity region, with respect to 
a fixed finite field, that inner bounds the corresponding network 
linear coding capacity region, show that it is a computable ratio- 
nal polytope, and provide exact algorithms and approximation 
heuristics. We show connections between computing these regions 
and a polytope reconstruction problem and some combinatorial 
optimization problems, such as the minimum cost directed Steiner 
tree problem. We provide an example to illustrate our results. 
The algorithms are not necessarily polynomial-time. 

I. Introduction 

In a seminal work in 2000, Ahlswede et al. H] introduced 
the network coding model to the problem of communicating 
information in networks. They showed that the extended 
capabilities of intermediate nodes to code on incoming packets 
give greater information throughput than in the traditional 
routing model and that the capacity of any multiple multicast 
network with a single source node is equal to the minimum 
of min-cuts between the source and receiver nodes. 

A few explicit outer bounds of capacity regions of networks 
exist; the max-flow/min-cut bounds is one such set of bounds, 
which were sufficient in the case of single-source multiple 
multicast networks. Harvey et al.fT] combined information 
theoretic and graph theoretic techniques to provide a com- 
putable outer bound on the network coding capacity regions. 
Yan et al. f3l gave an explicit outer bound that improves upon 
the max-flow/min-cut bounds and showed its connection to 
a minimum cost network coding problem. They used their 
results to compute the capacity region of a special class of 
3-layer networks. Thakor et al.[4| gave a new computable 
outer bound, based on characterizations of all functional 
dependencies in networks. Yan et al.|5l provided an exact 
characterization of the capacity region for general multi-source 
multi-sink networks by bounding the constrained region in the 
entropy space. However, they noted that explicitly evaluating 
the obtained capacity regions remains difficult in general. In a 
related work, Chan and Grant|j6l showed that even the explicit 
characterization of capacity region for single-source networks 



can be difficult since the computation of a capacity region 
involves the non-polyhedral set of entropy functions and that 
linear programming bounds do not suffice. 

The routing capacity region of networks is better understood 
via linear programming approaches. Cannons et al.|7| defined 
a notion of network routing capacity that is computable with a 
linear program and showed that every rational number in (0, 1] 
is the routing capacity of some solvable network. Yazdi et 
al.lSl and 191 extended a special case of Farkas Lemma called 
the "Japanese Theorem" to reduce an infinite set of linear 
constraints to a finite set in terms of minimal Steiner trees 
and applied the results to obtain the routing capacity region of 
undirected ring networks. In a subsequent work, Kakhbod and 
Yazdi[101 provided complexity results on the description size 
of the inequalities obtained and applied them to the undirected 
ring networks. 

In this paper, we study the network capacity region of 
networks along the lines of work by Cannons et al. Q- We 
define the network capacity region of networks analogously 
to the rate regions in information theory and show its notable 
properties when the finite field is fixed: closedness, bound- 
edness and convexity. In the case of routing, we prove that 
the network routing capacity region is a computable rational 
polytope and provide exact algorithms and approximation 
heuristics for computing the region. In the case of linear 
network coding, we define an auxiliary region, called the semi- 
network linear coding capacity region, which is a computable 
rational polytope that inner bounds the network linear coding 
capacity region, and provide exact algorithms and approxima- 
tion heuristics, with respect to a fixed finite field. The main 
idea is to reduce computation of the capacity regions to a 
polytope reconstruction problem and use linear programming 
techniques on associated combinatorial optimization problems. 
Our results generalize to directed networks with cycles and 
undirected networks. This present work partially addresses two 
problems proposed by Cannons et al. Q: whether there exists 
efficient algorithms for computing their notions of network 
routing capacity and network linear coding capacity. It follows 
from our work that there exist combinatorial approximation 
algorithms, not polynomial-time, for computing the network 
routing capacity and for computing a lower bound of the 
network linear coding capacity. In addition, we provide an 
example to illustrate our results. 



We note that algorithms and heuristics we present are not 
polynomial-time schemes. We do not distinguish the vertex 
and hyperplane descriptions of polytopes, but note that con- 
verting one description into another can be computationally 
expensive. 

The paper is organized as follows. In Section II, we give a 
fractional network coding model. In Section III, we define 
the network coding capacity region and prove its notable 
properties. In Section IV and V, we prove notable properties of 
network coding capacity regions, or inner bounds thereof, in 
the case of routing and linear coding and discuss algorithmic 
aspects of the regions. In Section VI, we give an example that 
illustrates our results. Finally, in Section VII, we discuss fur- 
ther extensions and conclude. Due to the page limit, we omit 
some details and refer to 111] for more detailed exposition. 

II. Fractional Network Coding Model 

We define a fractional network coding model. Most defini- 
tions are adapted from Cannons et al. [7|. We write vectors 
or points in a multi-dimensional space with a hat as in k and 
let ki (also k{i)) denote the i-th coordinate of the vector. We 
shall omit the hat when convenient. 

A capacitated network A/" is a finite, directed, acyclic 
multigraph given by a 7-tuple {v, e, /i, c, A, S, R) where i/ is 
a node set, e is an edge set, /i is a message set, c : e — > Z+ is 
an edge capacity function, A is an alphabet, S* : — > 2^ is a 
source mapping, and i? : 2^ is a receiver mapping. 

We define fractional edge function, fractional decoding 
function, and fractional message assignment with respect to 
a finite field F, where |F| > \A\, a source dimension vector 
k, and an edge dimension n. Let M = (v, e, /i, c, A, S, R) be 
a capacitated network and mi, . . . , m|p| be the messages. Let 
k = (fci, . . . , fc|^|) be a vector of positive integers and n be a 
positive integer. For each edge e = {x,y), a fractional edge 
function is a map /e : (F'^'i ) x • • • x (F'^^" ) x x 
... X (i^"c(e„+3)-j pnc(e)^ p^j. g^^j^ node X £ v and message 

rrij e R{x), ci fractional decoding function is a map fx,mj ■ 
(_pfcii)x- • •x(i^'='")x(i^"'=('="+i))x- • -x -> F''\ 

where m^^ , . . . , m^^ are a messages generated by x and 
Ba+i, ■ ■ ■ ,ea+j3 are f3 in-edges of x. We denote the collec- 
tions of fractional edge and fractional decoding functions by 
Te = {fe ■■ e € e} and Ta:,m = {fx.m ■ X e iy,m e R{x)}. 
Note ki is the source dimension for message rrii. 

A fractional message assignment is a collection of maps 
a — (ai, . . . , a|^|) where is a message assignment for 
rrii, Qi : nii F^* . A fractional network code on M is 
a 5-tuple {F,k,n,Ti,,Td) where F is a finite field, with 
\F\ > \A\, /c is a source dimension vector, n is an edge 
dimension, is a collection of fractional edge functions, and 
J^d is a collection of fractional decoding functions. We have 
different classes of fractional network codes corresponding 
to and being routing, Unear, and nonlinear functions 
over the field F: fractional routing/linear/nonlinear network 
codes. A fractional network code is a. fractional network code 
solution if for every fractional message assignment a, we 
have fx^^. {ui^ (m^ J, . . . , a^^ ("T-io ), s{ea+i), s{ea+p))^ 



aj{mj), for all x E v and nij E R{x), where rn^j^, . . . , m^^ 
are a messages generated by x and Ca+i, ■ ■ ■ , e^+p are /3 in- 
edges of X. If the above equation holds for a particular x E v 
and message m E R{x), then we say node x's demand m is 
satisfied. We have classes of fractional network code solutions: 
fractional routing/linear/nonlinear network code solutions. 

If {F, k, n, Je, Td) is a fractional network code solution for 
Af, source node x E ly sends a vector of ki symbols from F, 
representing alphabets in A, for each message rrii G S{x); 
each receiver node x E demands the original vector of ki 
symbols corresponding to rrii for each rrii E R{x); and each 
edge e carries a vector of c{e)n symbols. We refer to coor- 
dinates of the symbol vector corresponding to rrii as message 
rrii '-^ coordinates and coordinates of the symbol vector on edge 
e as edge e 's coordinates. Note that a coordinate of edge e can 
be active or inactive, depending on whether it actively carries 
a symbol in the fractional network code solution or not. A 
fractional network code solution is minimal if the set A of 
all active coordinates of edges is minimal, i.e., there exists no 
fractional network code solution with the same F, k and n 
and the set of active coordinates that is a strict subset of A. 

In this paper, we consider non-degenerate networks; for 
each demand of a message at a node, there is a path from 
a source node generating the message to the receiver node 
and no message is both generated and demanded by the same 
node. We shall abridge network code descriptions and solution 
concepts when convenient. In this work, a multicast network 
has exactly one message produced by a source node and 
demanded by more than one receiver nodes and multi-source 
multi-sink networks are the most general network instances. 

III. Network Capacity Regions 
A vector of nonnegative numbers E Q'^' 

is an achievable coding rate vector for M if there exists 
a fractional network code solution {F,k,n,J-e,J^d) where 
k = (fci, . . . , fc|^|). The network capacity region of is 
the closure of all achievable coding rate vectors in K'^L 
There are different classes of network capacity regions: the 
network routing capacity region, Cr, which is the closure of 
all achievable routing rate vectors; the network linear coding 
capacity region, Ci, which is the closure of all achievable linear 
coding rate vectors; and the network coding capacity region, 
C, which is the closure of all achievable coding rate vectors. 

Theorem 1: Let M a capacitated network and F a fixed 
finite field. The corresponding network capacity region C is a 
closed, bounded and convex set in 

Proof: (Closedness) By definition, C is a closed set. 
(Boundedness) By symmetry, it suffices to show that ^ is 
bounded in any achievable coding rate vector. Let n be the 
edge dimension, vi be the set of source nodes in v that 
generate message mi and 7 be the sum of capacities of out- 
edges of nodes in vi. Then, fci < as we cannot send 
more than 771 independent coordinates of message mi and 
expect receivers to recover all the information. Hence, ^ < 7. 
(Convexity) Let xq, xi G C and A E [0, 1]. It is straightforward 
to show that x — [1 — X)xq + Xxi e C by rate sharing; we 



use sequences of achievable coding rate vectors converging to 
xq and xi to produce a sequence of achievable coding rate 
vectors converging to x. ■ 
Similarly, the routing capacity region Cr and linear coding 
capacity region Ci are closed, bounded and convex regions 
when the finite field F is fixed. However, the rate-sharing 
argument seems to break down in the case of linear coding if 
the finite field is not fixed. 

IV. Network Routing Capacity Regions 

We show the network routing capacity region Cr is the 
image of a higher-dimensional rational polytope under an 
affine map and consider the computation of Cr as the polytope 
reconstruction problem with a ray oracle. Since multi-source 
multi-sink networks can be reduced to multiple multicast 
networks, it suffices to show the results with respect to the 
multiple multicast networks. 

A. Properties 

Theorem 2: The network routing capacity region C,. is a 
bounded rational polytope in m'^'' and is computable. 

Proof: (Polytope) It suffices to consider minimal frac- 
tional routing solutions, which consist of routing messages 
along Steiner trees, since any fractional routing solution can 
be reduced to a minimal one. Let % be the finite set of all 
Steiner trees rooted at the source node of message m,; and 
spanning all receiver nodes that demand rrii, and T be the 
union, T = Ti U . . . U 7[p|. Then, any minimal fractional 
routing solution satisfies the following constraints: 

^yg^ T(e) • a;(r) < c(e)n, Ve e e 

Y^TdTi ^(^) ^ VI < « < 

a; ' > 0, 

where x{T) is the number of times Steiner tree T is used 
in the solution and T(e) is an indicator that is 1 if T uses 
edge e, or otherwise. After scaling x{T) by n, any minimal 
fractional routing solution satisfies 

Y.TeTT{.e) ■ x{T) < c(e), Ve G e 
X > 0. 

As the coefficients are in Q, the above set of inequali- 
ties defines a bounded rational polytope Vr, with rational 
extreme points. The polytope is bounded, because edge ca- 
pacities are finite and no Steiner tree can be used to route 
for infinitely many times. Each minimal fractional routing 
solution reduces to a rational point inside the polytope Vr, 
and each rational point x inside Vr has a minimal fractional 
routing solution {F, k,n, J^^, J^d) that reduces to it, such that 
^ — X^TeT ■^(^) ^'^^ follows that the network routing 

capacity region Cr is the image of Vr under the affine map 

A ■■ {x{T))TeT ^ (ETeri^(^)'---'STe7]^,| ^(^))- 
follows that Cr is a bounded rational polytope. 

(Computability) We first compute the vertices , . . . , w/,, of 
polytope Vr by any vertex enumeration algorithm. Then we 
compute the images of the vertices of Vr under the affine 
map '4>r- The network routing capacity region is given by the 



vertices of the convex hull of points (wi ),..., V'r(i^h) in 

B. Algorithms 

We now provide exact algorithms and approximation heuris- 
tics for computing the network routing capacity region Cr- We 
already provided an exact algorithm in the proof of Theorem|2] 
Since the polytope Vr is defined in a high-dimensional space, 
the exact algorithm may not be efficient in practice. Using 
results in |12| and |13| on polytope reconstruction problems, 
we derive different kinds of algorithms that might be more 
efficient. 

In our polytope reconstruction problem, we compute the 
facet description of a given polytope by making calls to the 
ray oracle Onay which, given a ray, returns the intersection 
point on the polytope surface and the ray. We reduce the 
computation of Cr to the polytope reconstruction problem by 
1) reflecting Cr around the origin to get a symmetric polytope 
Q that contains the origin in its interior and 2) solving the 
linear programs similar to the one in Cannons et al. |7| to 
implement the ray oracle Ouay To reflect Cr, we map all 
calls to the ray oracle to equivalent calls with rays defined in 
m'[''. We use the algorithm outlined in Section 5 of Gritzmann 
et al. (13] to compute all the facets of the resulting polytope 
Q and, therefore, Cr- The main idea of the algorithm is to first 
find a polytope Q' that contains Q and whose facet-defining 
hyperplanes are a subset of those for Q (Theorem 5.3 in fT3\), 
and then successively add more facet-defining hyperplanes of 
Q to Q' by using Ouay- By Theorem 5.5 in Gritzmann et 
al. IIT3I and the symmetries around the origin, we need at most 
MCr) + (ImI - + (5|m| -4)/|^|_i(C) calls to the 

ray oracle where fi{Cr) denotes the number of i-dimensional 
faces of Cr that do not contain the origin (the 0-th dimensional 
faces being the points). 

If we use an exact algorithm for the ray oracle Ojiay, we 
get an exact hyperplane description of the network routing 
capacity region. If instead we use an approximation algorithm 
that computes some point r such that the actual intersection 
point lies between r and Ar, then we obtain an approxima- 
tion heuristic that computes a set of points r such that the 
boundary of Cr lies between points r and Ar- We note that 
an approximation algorithm for Opiay does not necessary yield 
an approximation algorithm forCr, where an A-approximation 
of Cr would be a polytope V such that V C Cr C AV- 
While an approximation algorithm for the oracle Ojiay does 
not necessarily lead to a polytope description of Cr, it might 
be faster and more efficient than exact algorithms and, hence, 
more applicable to compute a quick "sketch" of Cr- For 
instance, we can find approximate intersection points on a 
sufficiently large number of rays evenly spread apart. 

C. Implementations of Oracle Ojiay 

Given the hyperplane description of the polytope Vr and a 
ray with a rational slope, x = qt,t > 0, we want to compute 
the rational intersection point of the ray and the boundary of 



C,.. Note that the intersection point is exactly 
^max is the optimal value to the linear program 



where 



max A 
s.t. X^TeT 

X, X 



T{e)-x{T) < c{e), Ve € 6 

> Xqi, \fi 

> 0. 



We can use any linear programming algorithm, such as 
the ellipsoid algorithm, to solve the linear program exactly 
and, thus, obtain an exact oracle Oriay For the approximate 
ray oracles, we can employ algorithms that solve the linear 
program within a provable approximation guarantee. Note 
that as the network routing capacity region Cr is a rational 
polytope, it suffices to consider rays with a rational slope. 

Using the results for the multicommodity flow and related 
problems by Garg and Konemann lfT4l . we can derive a 
combinatorial approximation algorithm that solves the linear 
program ( |IV. It computes a point f such that XmaxQ is 
on the line segment between f and (1 + ijj)Af for some 
numbers w > and A > 1. The main idea is to view 
solving ( |IV. 1 | i as concurrently packing Steiner trees according 
to the ratio defined by q and use the results for the minimum 
cost directed Steiner tree problem: given an acyclic directed 
multigraph Q = {v, e), a length function / : e — > M+, a 
source node s and receiver nodes ni, . . . ,nk, find a minimum 
cost (defined by Xeec' ^(^)) subset of edges e' such that 
there is a directed path from s to each Ui in e'. We assume 
we have an oracle Onsteiner that solves the minimum cost 
directed Steiner tree problem, which is well-known to be NP- 
hard, within an approximation guarantee A. Then we have the 
following theorem: 

Theorem 3: There exists an (1 + a;)^-approximation algo- 
rithm for the linear program ( |IV. l| i in time 0{oj^^ log Aj/ij 
+ |e|)^log|e| • Tosteiner), whcrc Tosteiner is flie time re- 
quired to solve the minimum cost directed Steiner tree problem 
with oracle Od steiner within an approximation guarantee A. 

Note that computing the network routing capacity region is 
not the same as the multicommodity flow, as in our problem, 
a message can be demanded by multiple nodes and can 
be duplicated at intermediate nodes. There are algorithms 
for Ojjsteiner', the approximation using shortest paths yields 
A = 0{\u\) in time 0((|!^p + |i^||e|) log and Charikar et 
al. ifTSl gives a set of algorithms with A = i{i — 
and time 0(|;/p*) for any integer i > 1. We can also 
use any brute-force algorithm that yields A = 1. Note that 
further improvements in the minimum cost directed Steiner 
tree problem translates to improvements to our approximation 
algorithm by above theorem. 

V. Network Linear Coding Capacity Regions 

We define a computable polytope C;, which we call the 
semi-network linear coding capacity region, that is contained 
in the network linear coding capacity region C;. It is unknown 
how good of an approximation the polytope C[ is to C;. 
Unlike in the last section, the finite field is important in 
the computation of C[ because of linearity of functions and 



termination of algorithms. In this section, we assume that a 
network M is a. multi-source multi-sink network and that the 
finite field F is fixed. 

A. Definitions 

The weight vectors associated with Af are vectors w in 
{0, such that there exists a scalar-linear network code 
solution when only messages rrii with Wi — 1 are considered 
and all the edges are assumed to be of unit capacity. We refer 
to the scalar-linear network code solutions corresponding to 
these weight vectors as partial scalar-linear network code 
solutions. A fractional network code is a simple fractional 
linear network code solution if the fractional network code 
is linear over the finite field F and can be written as an 
aggregate of partial scalar-linear solutions (when considered 
with unit edge capacities). We define the semi-network linear 
coding capacity region C[ of M as the closure of all coding 
rate vectors achievable by simple fractional linear network 
code solutions. Clearly, the network linear coding capacity 
region C/ contains the semi-network linear coding capacity 
region C[. Note that scalar-linear network code solution and 
simple fractional linear network code solution are two different 
solution concepts. 

B. Results 

By the same line of reasoning as in the case of routing, we 
get the following theorems. 

Theorem 4: Assume a finite field F is given. The semi- 
network linear coding capacity region €[, with respect to F, 
is a bounded rational polytope in m'|'' and is computable. 

Instead of the minimum cost directed Steiner tree problem, 
we have two associated problems: minimum cost scalar-linear 
network code problem where given a network with unit edge 
capacities, a finite field F, and a length function / : e — > M+, 
we want to compute the minimum cost (defined to be the 
sum of lengths of the edges used) scalar-linear network code 
solution; and fractional covering with box constraints problem 
which is to solve 

min c{j)x{j) 

s.t. EjMhjMj) > b{i), yi<i<n 

x{j) < u{j), Vl<j<m 

X > 0, 

where A is an n x m nonnegative integer matrix, b a non- 
negative vector, c a positive vector, and u a nonnegative 
integer vector Assuming we have two oracles, OsLinear and 
Tpcover, for thcsc problems, we get: 

Theorem 5: There exists an (1 + ti;)i?-approximation algo- 
rithm for oracle O^ay, for C;, with time 0{uj^'^ {log A\fi\ + 

\€\)B log |e| • {TpCover + k' Ts Linear)) , wherC TpCover is the 

time required to solve the fractional covering problem by 
Opcover within an approximation guarantee B, Ts Linear is 
the time required to solve the minimum cost scalar-linear 
network code problem exactly by Os Linear, and k' is the total 
number of weight vectors associated with M. 
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Fig. 1. Network and its network routing capacity region and semi- 
network linear capacity region, eacli computed witli an exact Onay and an 
approximate Onay 

For Opcover, wc caii use any linear programming al- 
gorithm or the combinatorial approximation algorithm by 
Fleischer ||T6| . For Os Linear, there is a dynamic programming 
algorithm in [11|. Without the minimum cost condition, the 
scalar-linear network code problem reduces to the decidability 
problem of determining whether or not a network has a scalar- 
linear solution, which is NP-hard WH . Without the fixed finite 
field F, the decidability of the problem is unknown; this 
justifies in part our assumption of the fixed finite field F. 

VI. An Example 

For network M, we computed the network routing capacity 
region (two inner curves) and the semi-network linear coding 
capacity region (two outer curves) with both exact and ap- 
proximate ray oracles in Figure [T] Network M has two source 
nodes at the top and two receiver nodes at the bottom. For 
exact ray oracles, we hard-coded the corresponding linear pro- 
grams and used a linear program solver, linprog, in MATLAB. 
To obtain approximate intersections points, we simply used 
an approximate oracle Onay in place of the linear program 
solver in an implementation of a 2D-polytope reconstruction 
algorithm. We note that the approximate oracle Opiay worked 
well in this example and led to its successful termination, but 
this may not hold for arbitrary networks in general. 

VII. Further Discussion 

Our results have a few extensions: We can design mem- 
bership testing algorithms that given a rate vector, determines 
whether or not there exists a fractional network code solution 
that achieves the rates from algorithms we have provided, for 
the network routing capacity region and semi-network linear 
coding capacity region. We can also generalize the results 
to directed networks with cycles and undirected networks 
straightforwardly. 

Note that the network routing capacity defined by Cannons 
et aLQ corresponds to a point on the boundary of poly tope C^; 
it is exactly the intersection point between the (outer) boundary 
of Cr and the ray x — (1, . . . , l)t,t > 0. Hence, our work 
partially addresses two problems proposed by Cannons et al.: 
whether there exists efficient algorithms for computing their 



notions of network routing capacity and network linear coding 
capacity. It follows from our work that there exist combinato- 
rial approximation algorithms, albeit not polynomial-time, for 
computing the network routing capacity and for computing a 
lower bound of the network linear coding capacity. 
For details omitted in this paper, we refer to 1, 1 1 J . 
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